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Abstract 

We report a measurement of the inclusive electron energy spectrum for semileptonic decays of 
B mesons in a 140 fb~^ data sample collected near the T(4S') resonance, with the Belle detector at 
the KEKB asymmetric energy e'^e~ collider. We determine the first and second moments of the 
spectrum for threshold values of the electron energy between 0.6 and 1.5 GeV. 

PACS numbers: 12.15.Hh, ll.30.Er, 13.25.Hw 
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INTRODUCTION 



The Cabibbo-Kobayashi-Maskawa matrix element \Vcb\ can be extracted from the inclu- 
sive branching fraction for semileptonic b hadron decays B{B Xcii')^,^. Several studies 
have shown that the spectator model decay rate is the leading term in a well-defined ex- 
pansion controlled by the parameter AqcTi/'^b- Non-perturbative corrections to this leading 
approximation arise only at order and above. The key issue in this approach is the 

ability to separate non-perturbative corrections, which can be expressed as a series in powers 
of l/rrih, and perturbative corrections, expressed in powers of a^. 

The coefficients of the l/rrih power terms are expectation values of operators that in- 
clude non-perturbative physics. Two different expansions exist, reflecting a difference in 
the approach used to handle the energy scale fj, that separates long-distance from short- 
distance physics. In this note, we use the short distant mass expansion that defines the 
non-perturbative operators using a mass scale fi ^ 1 GeV 0. 

The shape of the lepton spectrum provides constraints on the heavy quark expansion 
based on local Operator Product Expansion (OPE), which calculates properties of the 
B — *■ Xcii' transitions. So far, measurements of the electron energy distribution have been 
made by the DELPHI, CLEO, and BABAR collaborations UMM- 

In this note we report 

a measurement of the first and second moment of the electron energy spectrum with a min- 
imum electron momentum cut ranging between 0.6 and 1.5 GeV in the B meson rest frame. 
We measure, independently, the electron energy moments for the semileptonic decays of 
the B~^ and the B^ mesons 0. If quark hadron duality applies in 5 — > Xcii' decays then 
the moments from the B~^ and B^ should agree. Up to now such a measurement has been 
impossible to perform. However, in this measurement we find good agreement between the 
moments. 

The data used in this analysis was collected with the Belle detector at the KEKB 
asymmetric energy e~^e~ collider. The Belle |2] detector is a large-solid-angle magnetic spec- 
trometer that consists of a three-layer silicon vertex detector (SVD), a 50-layer central drift 
chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), a barrel-like ar- 
rangement of time-of-fiight scintillation counters (TOP), and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a super-conducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron fiux-return located outside of the coil is instrumented 
to detect mesons and to identify muons (KLM). 

Events are selected by fully reconstructing one of the B mesons, produced in pairs from 
T(4S'), in several hadronic decay modes. Prompt semileptonic decays {b xiu) of the 
non-tag side B mesons are separated from cascade charm decays {b ^ c ^ yC-^)-, based on 
the correlation between the fiavour of the tagged B and the lepton charge. 

The present results are based on a 140 fb~^ data sample collected at the T(4S') resonance, 
which contains 1.5 x 10*^ BB pairs. An additional 15fb~^ data sample taken at 60 MeV 
below the T(4S') resonance (off— resonance) is used to perform background subtraction from 
the e^e~ qq process. 

We use a fully simulated generic Monte Carlo sample generated with the qq98 event 
generator This sample is equivalent to about three times the beam luminosity of the 
real data. Simulated events are required to satisfy tight hadron selection criteria. 
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FULL RECONSTRUCTION SAMPLE 



The hadronic decay of the tag-side B meson is fully reconstructed in the decay modes 
B D^*h+,D^*^p+,D^*1ai, yielding a high purity B meson sample. The following sub- 
decay modes are implemented: 

. D*+ ^ D\+,D+7r^, 

• Kit, Ktttt^, KinnT, Kstttt, Kstc^ and 

• — > Ktttt, Kstt. 

If multiple candidates are found in one event, we select the best candidate based on AE 
and other variables. There are two conditions that affect the best candidate selection. The 
candidate may be reconstructed from at least one particle from the other side B, or the final 
state particles may be incorrectly reconstructed. 

ELECTRON SELECTION 

We identify electrons produced in semileptonic B decays on the "non-tag" side. Low- 
momentum particles spiral in the CDC detector. If the particles pass close to the interaction 
point more than once, they can lead to multiple reconstructed tracks. Those tracks are 
removed. The remaining tracks are required to satisfy tight restrictions on the impact 
parameters. In addition, we require the tracks to be in the barrel region of the detector, 
corresponding to an angular acceptance of 35° < ^lab < 125°, where 6'iab is the polar angle 
of the track relative to the z axis (opposite the positron beam line). 

Tracks that pass the above selection criteria and are not used in the reconstruction of 
the tag-side B meson are considered as electron candidates. Electron identification is based 
on a combination of specific ionisation dE/dx measurements in the CDC, the response of 
the ACC, the shower shape in the ECL, and the ratio of energy deposit in the ECL to 
the momentum measured by the tracking system (E/p) |^. The electron identification 
efficiencies include the detector acceptance as well as the selection efficiency. The errors for 
the tracking and electron detection efficiencies are ±1% and ±2%, respectively. Electron 
identification efficiency is determined by analysing a data sample where simulated single 
electron tracks are overlaid on hadronic events taken from real experimental data. The 
electron momentum spectrum is corrected by using this momentum dependent electron 
detection efficiency. 

The momentum of electron candidates is measured in the B meson rest frame (p**^^-*), 
exploiting the knowledge on the momentum of the fully reconstructed B. To ensure a good 
electron identification performance we require p^^^^ > 0.6 GeV/c. 

Backgrounds from J/ip, V^(25'), photon conversions in the material of the detector, and 
Dalitz decays from 7r° and rj mesons are reduced by imposing veto cuts. We look for addi- 
tional electrons in the event, which satisfy looser identification criteria and Pe > 0.3 GeV/c. 
We then calculate invariant masses {rriee) for each electron candidate when combined with 
an opposite-charge electron that passes the loose selection. In the Dalitz case we require 
an additional photon (mge^). We reject the electron if m^e lies within the nominal J/ip and 
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ij^i^S) mass regions, or mee^y lies within the nominal 7r° and rj mass regions. To remove 
photon conversions we reject candidates by placing cuts on the vertex of the decay and rriee 
less than 100 MeV/c^ 

NUMBER OF TAGGED EVENTS 

We select events containing a fully reconstructed B candidate on one side and an electron 
candidate on the other side. 

For each selected event, we calculate the beam-constrained mass Mbc and AE: 

Mbc = v/(^LanJ' - iPhV, AE = E*^ - E;,,^, (1) 

where -Ebeam' beam energy, the B 3-momentum and the B energy in 

the centre of mass frame, respectively. Events with Mbc > 5.27 GeV/c^ and —0.06 GeV < 
AE < 0.08 GeV are considered to be signal candidates. 

Each electron is assigned to the "right sign" sample if the electron has the opposite 
flavour to the fully reconstructed B, and to the "wrong sign" sample if the electron has 
the same flavour as the fully reconstructed B. In events without B^B'^ mixing, the primary 
electrons belong to the "right sign" sample while the secondary electrons contribute to the 
"wrong sign" sample. Figure ^ shows the beam constrained mass distribution for B^ and 
B~^ right sign and B^ wrong sign electrons, with scaled off resonance background overlaid. 
The number of events in the signal region, is 3367 ± 58 (1659 ± 40) and 6831 ± 83 for the 
B^ right sign (wrong sign) and B^ right sign candidates, respectively. 

BACKGROUND SUBTRACTION 

The selected electron energy spectrum is contaminated by background processes, which 
should be evaluated and subtracted from the distribution before the extraction of the mo- 
ments. The residual background is from: 

• continuum background, 

• full-reconstruction background, 

• background from secondary decays, 

• J/ip,ip{2S), Dalitz decays and photon conversions where one electron of the pair has 
escaped detection and 

• fake electrons. 

The electron momentum spectrum of the continuum background is derived from 
off— resonance data and is normalised using the off— to on— resonance luminosity ratio. 
To account for the low statistics in the off— resonance data we fit an exponential to the 
electron energy distribution. We fit only the B~^ candidates then apply the same continuum 
cut to the -B°, scaled to the electron yield of the wrong sign and right sign B^ samples. 
There are insufficient statistics to fit the off— resonance B^ samples. 

The electron momentum spectrum of the combinatorial background is derived from the 
generic BB Monte Carlo events where either reconstruction or flavour assignment of the 
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FIG. 1: The beam constrained mass after electron selection cuts and AE cuts for the right sign 
electron sample (a), right sign electron sample (b) and the B^ wrong sign electron sample (c). 



tagged B meson is not carried out correctly. The BB Monte Carlo events are normalised 
to the on resonance data Mbc side band (Mbc < 5.25 GeV/c^) after continuum background 
subtraction in this region. 

We also correct for cases where the fully reconstructed B is correctly tagged, but the 
electron does not directly originate from a B decay or is a hadron misidentified as an 
electron. These backgrounds are irreducible; to estimate the magnitude of the contributions 
we normalise the BB Monte Carlo to the electron yield after continuum and combinatorial 
background subtraction. 

Secondary electrons arising from B D^*^ e, B ^ r ^ e, and B ^ D^*^ ^ r ^ e 
decays are simulated using the latest published branching fractions [l^ . 
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Contributions from J/V', V'(25') decays, photon conversions, and Dalitz decays are insignif- 
icant after the electron selection cuts. The surviving backgrounds, where one electron of 
the pair has escaped detection, are estimated by the Monte Carlo simulation and subtracted 
along with the major secondary backgrounds. 

The Monte Carlo yield due to fakes is corrected for the difference in the fake rate for 
data and Monte Carlo measured with samples of tt^tc^ decays. The magnitude of 

the correction is assigned as the systematic uncertainty on the fake yield. 

Figure El shows the raw electron momentum spectrum with all background contributions 
overlaid. Table H] summarises the number of detected electrons and the contributions from 
these backgrounds. 



MIXING CORRECTIONS 

To reduce the contamination from cascade semileptonic charm decays, we determine 
whether the electron's parent and the fully reconstructed B have the same or opposite 
flavour. For events tagged as we require the right sign tag. For events tagged as B^, 
we introduce a dependence on B^B^ mixing which is taken into account by solving the 
equations; Alright = A^p(l - Xd) + N^Xd, A^wrong = NpXd + A^c(l - Xd) where Alright and A^wrong 
are the number of electrons with the right and wrong signs, respectively. A^p and Nc are 
the number of electrons from the prompt and cascade semileptonic decays, respectively and 
Xd = 0.186 ± 0.004 fl^ is the mixing probability. 



CORRECTIONS TO THE ELECTRON ENERGY SPECTRUM 

The electron energy spectrum is generated via Monte Carlo simulation of B ^ XcCu 
decays using the qq98 event generator |lO|. The spectrum from B XcCu is modelled 
using four components: = D fISGW2 D* (HQET higher resonance charm 



meson states D**(ISGW2) and non-resonant D^*^tc (Goity and Roberts [13|). To account for 
the most recent theoretical and experimental results, we re-weight the D and D* components 
in to the spectra generated with current world average form factors [l^ . 

The available Monte Carlo sample does not incorporate 0{a) QED corrections. We 
correct the background subtracted electron spectrum with energy-dependent weights using 
the PHOTOS O package. 

Electrons that come from the h ^ u transition are subtracted from the unfolded electron 
energy spectrum, defined later. We model the electron energy spectrum from the B — >■ XJ.i' 
using the De Fazio and Neubert prescription . The 6-quark motion parameters are 
derived in Ref. We scale according to the B Xulv branching fraction in Ref. Ei. 



The electron momentum spectra, after all corrections have been applied, are shown in 
Figure 01 



UNFOLDING THE ELECTRON ENERGY SPECTRUM 

To measure the first and second electron moments we need to determine the true elec- 
tron energy spectrum. The background subtracted energy spectrum is distorted by various 
detector effects. The true electron energy spectrum is extracted by performing an unfolding 
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FIG. 2: Breakdown of the backgrounds in the electron momentum spectra for right sign (a), 
right sign (b) and wrong sign (c) electrons. 



procedure based on the Singular Value Decomposition (SVD) algorithm of A.Hocker and V. 
Kartvelishvili [l8| . 

Figure m shows the correlation between the reconstructed and generated electron energy. 
Due to the good resolution of the Belle detector, the reconstructed electron energy is highly 
correlated with the generated one. The off— diagonal bins are filled due to bremstrahlung 
occurring in the material in front of the calorimeter. It is important to note, for the reliability 
of the unfolding, that the phase space covered by the generated energy spectrum is very 
similar to the reconstructed spectrum. 

We divide the generated energy in slices of 0.025 GeV. For each slice we produce a 
histogram of the corresponding reconstructed energy spectrum, normalised to one. The 
histogram has a bin width of 0.025 GeV. Bin and slicing widths have been chosen to have 
less than 40% correlation (cross-talk) between bins and/or slices. These histograms are then 
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TABLE I: Number of Electrons 



B candidate 


right-sign 


B^ wrong-sign 


5+ right-sign 


On Resonance Data 


3367± 58 


1659 ± 40 


6831 ± 83 


Scaled Off Resonance 
Combinatorial Background 


54 ± 7 
420 ± 26 


30 ± 6 
339 ± 29 


106 ± 10 
1073 ± 86 


Bn Background 
Secondary 
Hadron Fakes 


479 ± 21 
456 ± 20 
23 ± 4 


856 ± 29 
847 ± 29 
9 ± 2 


688 ± 26 
634 ± 23 
54 ± 5 


Background Subtracted 
After Efficiency Correction 
After Mixing Correction 


2414 ± 67 
4066 ± 115 
5054 ± 145 


434 ± 58 
727 ± 100 

N/A 


4964 ± 122 
8371 ± 209 
N/A 
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Pg*(B) (GeV/c) Pg*(B) (GeV/c) 

FIG. 3: The electron momentum spectrum in the B rest frame after background subtraction, 
mixing and efficiency correction for B^ (a) and B^ (b). 

used to generate a response matrix. Next, we use the SVD algorithm to perform singular 
value (eigenvalue) analysis of the response matrix. Figure El compares the unfolded Monte 
Carlo electron energy spectrum compared to the true spectrum. 

We weight the Monte Carlo spectra with the QED radiative effects in order to calculate 
the correct unfolding matrix. When we unfold the data we remove the QED radiative effects, 
as the OPE does not have an 0{a) QED correction. The unfolded electron energy spectrum 
for data is shown in Figure IHl 

MOMENTS OF THE ELECTRON ENERGY SPECTRUM 

The first moment of the electron energy spectrum is defined to be (E'cut) = (-^e) E^>E^^t 
and the second M|(i?cut) = ((-^e — Ml{Ecnt)Y) Ee.>E^^,t- We measure the first and second 
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FIG. 4: Correlation between reconstructed and generated electron energy (i.e. reconstructed versus 
generated electron energy spectrum) for (a) and (b). The high population in the diagonal 
bins demonstrates the high resolution in the unfolding procedure. 
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FIG. 5: Generated and unfolded distributions for Monte Carlo simulated events (b). The 
reconstructed electron spectrum after the detector simulation is show in (a). The continuous line 
in (b) is the true (generator level) Eg distribution. The unfolded distribution is overlaid (data 
points). 



moments with five electron energy threshold cuts (-Ecut = 0.6, 0.8, 1.0, 1.2 and 1.5 GeV) 
for the and mesons. Table IIVI gives the final measurements of the moments, with 
radiative corrections applied. 

Figure [7| shows the measured moments for the and B^ data samples with the the- 
oretical bounds as a function of the threshold electron energy. The theoretical bounds in 
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FIG. 6: Unfolded Ee distribution corrected for QED radiative effects for (a) and (b). The 
errors shown are statistical. 



Figure [7| are taken from Ref. [3|, in the kinetic mass scheme. 



SYSTEMATIC UNCERTAINTIES 

The systematic uncertainty in the moments stems from event selection, electron identifi- 
cation, background estimation, model dependence, and mixing. 

The uncertainty due to mis-tagging in the and samples is derived from the mag- 
nitude of the combinatorial background and is estimated to be ±2%. Uncertainty in the 
measured luminosity provides an uncertainty on the continuum electron yield of ±1%. 

The effect of veto cuts has been calculated as the quadratic sum of averaged variations in 
the moments when selection criteria are loosened and tightened about their chosen positions. 

Subtraction of b —>■ u decays acquires an uncertainty associated with the renormalisation 
to the full inclusive spectra. The renormalisation factor is varied by ±50%. 

The uncertainty in secondary [B ^ D ^ e) decays derives from uncertainty in the value 
of the branching fractions. Model dependence is estimated from the observed change in the 
moments when the B — >■ Dev and B — > D*ev decay shape parameters are varied. 

The uncertainty in the mixing probability Xd is based on the current quoted error in Ref. 

a. 

In addition, we calculate systematic uncertainty associated with the magnitude of the 
hadron fake contribution, electron tracking efficiency, and electron detection efficiency. 

The total systematic error is obtained by adding each contribution in quadrature. It is 
important to note that the systematic errors are limited by Monte Carlo statistics, most 
significantly in the B^ sample. The contributions to the systematic error are summarised in 
table nil and table Hill 
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TABLE II: Breakdown of the systematic errors for the moments, Mi, M2 for — > X^eiy in the 



B meson rest frame for 3 values of E'cut 





Ml 


Ml 


Ml 


M2 


M2 


M2 




[MeV] 


[MeV] 


[MeV] 


[lO^MeV^j 


[lO^MeV^] 


[lO^MeV^] 


^cut[GeV] 


0.6 


1.0 


1.5 


0.6 


1.0 


1.5 


conversions and Dahtz decays 


0.02 


0.12 


0.04 


0.07 


0.00 


0.00 


e from J/tJj or ip{2S) 


0.91 


0.67 


0.63 


0.32 


0.15 


0.06 


B I?;:; ^ e 


1.89 


0.43 


0.02 


0.65 


0.11 


0.00 




0.02 


0.01 


0.00 


0.00 


0.00 


0.00 


^ u subtraction 


0.17 


0.18 


0.19 


0.09 


0.08 


0.06 


continuum background 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


combinatorial background 


1.58 


0.57 


0.01 


0.42 


0.10 


0.00 


hadron fakes 


1.24 


0.20 


0.01 


0.47 


0.06 


0.00 


electron detection 


0.01 


0.01 


0.01 


0.00 


0.01 


0.00 


tracking efficiency 


0.00 


0.0 


0.0 


0.0 


0.0 


0.0 


total systematics 


3.59 


1.13 


0.67 


1.21 


0.27 


0.09 



TABLE III: Breakdown of the systematic errors for the moments, Mi, M2 for B^ — > Xccv in the 



B mcHoii rc-st fraiiic^ for 3 x'ahu's of E^■^^l 





Ml 


Ml 


Ml 


M2 


M2 


M2 




[MeV] 


[MeV] 


[MeV] 


[lO^MeV^j 


[lO^MeV^] 


[lO^MeV^] 


^cut[GeV] 


0.6 


1.0 


1.5 


0.6 


1.0 


1.5 


conversions and Dalitz decays 


0.02 


0.26 


0.17 


0.21 


0.06 


0.01 


e from J/V' or ip{2S) 


0.40 


0.25 


0.51 


0.45 


0.20 


0.01 


B ^ D^;l ^ e 


2.38 


0.77 


0.05 


0.71 


0.13 


0.00 


X 


0.02 


0.00 


0.01 


0.01 


0.01 


0.00 


£)(*) 


0.18 


0.03 


0.00 


0.07 


0.01 


0.00 


b ^ u subtraction 


0.38 


0.37 


0.37 


0.18 


0.16 


0.12 


continuum background 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


combinatorial background 


0.85 


0.31 


0.14 


0.20 


0.00 


0.02 


hadron fakes 


0.93 


0.17 


0.01 


0.34 


0.03 


0.00 


electron detection 


0.12 


0.01 


0.01 


0.06 


0.01 


0.00 


tracking efficiency 


0.06 


0.00 


0.01 


0.03 


0.00 


0.01 


total systematics 


2.76 


0.99 


0.67 


0.97 


0.30 


0.12 



SUMMARY 

We report a measurement of the electron energy spectrum of the inclusive decay B — > 
XcCu and its first and second moments for threshold energies from 0.6 GeV to 1.5 GeV. This 
set of moments, combined with the measurements of the semileptonic branching fraction and 
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TABLE IV: Measured moments, Mi, M2 for B — > XcCv in the B meson rest frame for five cutoff 
energies -Bcut- Tlie first error is the statistical, and the second error is the systematic. The moments 
are corrected for QED radiative effects (using the PHOTOS package). 



i?cut[GeV] 


Ml [MeV] 


Ml [MeV] 


M2[103MeV2] 


M2[103MeV2] 




B+ 


B^ 






0.6 


1432.1 ± 4.3 ± 3.6 


1444.9 ± 5.5 ± 2.8 


150.1 ± 1.8 ± 1.2 


144.0 ± 2.1 ± 1.0 


0.8 


1487.2 ± 3.9 ± 2.2 


1488.0 ± 5.1 ± 1.8 


118.4 ± 1.4 ± 0.7 


119.0 ± 1.8 ± 0.6 


1.0 


1554.1 ± 3.6 ± 1.1 


1551.5 ± 4.7 ± 1.0 


88.1 ± 1.1 ± 0.3 


90.7 ± 1.4 ± 0.3 


1.2 


1631.7 ± 3.3 ± 0.7 


1632.6 ± 4.3 ± 0.8 


61.7 ± 0.8 ± 0.1 


64.1 ± 1.1 ± 0.2 


1.5 


1774.8 ± 2.8 ± 0.7 


1778.2 ± 3.8 ± 0.7 


30.6 ± 0.5 ± 0.1 


32.3 ± 0.7 ± 0.1 




Ecu,(GeV) E,,,(GeV) 



FIG. 7: The first and second electron energy moments. Mi (a) and M2 (b), as a function of cutoff 
energy E'cut with the bounds on the theory predictions shown as curved lines. The errors shown 
are statistical and systematic. 

the moments of the hadronic mass distribution, will be used for the determination of the 
HQE parameters and of \Vcb\- 
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